arXiv:1505.05857vl [physics.optics] 21 May 2015 


Observation of molecular dipole excitations by attosecond self-streaking 

Georg Wachter^0 Stefan Nagele^, Shunsuke A. Sato^, Renate Pazourek^, Michael Wais^, 
Christoph Lemelh, Xiao-Min Tong^’^, Kazuhiro Yabana^’^, and Joachim Burgdorfer^’® 

^Institute for Theoretical Physics, Vienna University of Technology, IO 4 O Vienna, Austria, EU 
^Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, Japan 
^Center for Computational Sciences, University of Tsukuba, Tsukuba 305-8577, Japan 
‘^Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Lousiana, USA and 
^Institute of Nuclear Research of the Hungarian Academy of Sciences (ATOMKI), Debrecen H-fOOl, Hungary, EU 

(Dated: May 22, 2015) 

We propose a protocol to probe the ultrafast evolution and dephasing of coherent electronic 
excitation in molecules in the time domain by the intrinsic streaking field generated by the molecule 
itself. Coherent electronic motion in the endohedral fullerene Ne@C6o is initiated by a moderately 
intense femtosecond UV-VIS pulse leading to coherent oscillations of the molecular dipole moment 
that persist after the end of the laser pulse. The resulting time-dependent molecular near-field 
is probed through the momentum modulation of photoemission from the central neon atom by a 
time-delayed attosecond XUV pulse. Our ab-initio time-dependent density functional theory and 
classical trajectory simulations predict that this self-streaking signal accurately traces the molecular 
dipole oscillations in real time. We discuss the underlying processes and give an analytical model 
that captures the essence of our ab-initio simulations. 


In recent years, the availability of wave-form controlled 
few-cycle intense laser pulses has afforded novel oppor¬ 
tunities to investigate the ultrafast and non-linear op¬ 
tical response of matter. Studies of electronic motion 
in the time domain have focused on rare-gas atoms [l| 
and molecules Si and, more recently, on nanostruc¬ 
tures, surfaces, and bulk matter Si. With the advent 
of attosecond XUV pulses synchronized to intense optical 
pulses, attosecond streaking Si has yielded insight into 
the real-time motion of electrons during photoionization 
i . Attosecond surface streaking S[i3 has found photo¬ 
electron spectra to be highly sensitive to the atomic-scale 
surface polarization and near-held distribution. These re¬ 
sults nurture the hope that attosecond streaking can also 
extract information on induced dipoles and near-fields in 
nano-scale systems bridging the gap between single atoms 
and extended matter. One chal leng e for such proposed 
“nano-plasmonic streaking” ll|, [l2| is the smearing out 
of the streaking signal over the nanometric target includ¬ 
ing emission from deeper layers. 

In this work, we analyze a scenario that allows to probe 
molecular coherent dipole excitation with sub-nanometer 
spatial resolution and attosecond time resolution by self¬ 
streaking. As prototypical case we study the endohedral 
fullerene Ne@C6o[l3- The Cgo cage serves as a nano¬ 
sized system with strong near-fields for which an ab-initio 
treatment is still feasible. The central neon atom serves 
as a distinct atomic-scale localized source of electrons to 
be ionized by an attosecond XUV pulse that probe the 
motion of the cage electrons. While transient molecular 
dipole moments during the IR laser pulse have recently 
been studied by Neidel et al. 1^, we employ here UV- 
VIS radiation to access low-lying electronic excitations 
of the molecule. The resulting coherent wavepacket mo¬ 


tion leads to an oscillating time-dependent dipole mo¬ 
ment sustained after the conclusion of the UV-VIS laser 
pulse. This time-dependent dipole moment causes an ap¬ 
preciable intrinsic streaking field for photoelectrons emit¬ 
ted from the central atom by a delayed attosecond XUV 
pulse. The streaking signal is found to trace the instan¬ 
taneous molecular dipole moment after the laser pulse, 
enabling the study of electronic excitations, their evolu¬ 
tion, dephasing, and decay with sub-femtosecond time 
resolution. 


In standard attosecond streaking the attosec¬ 

ond XUV pulse acts as “pump” which ionizes an elec¬ 
tron whose emission time is probed by the simultane¬ 
ously present near-IR probe pulse. In the present sce¬ 
nario, their roles are reversed: the UV-VIS pulse acts 
as a pump to excite the molecule, and the XUV pulse 
probes the molecular near-field associated with this ex¬ 
citation by a momentum shift of the outgoing electron 
(“self-streaking”) somewhat resembling attosecond tran¬ 
sient absorption [IMl- It can be viewed as a variant 
of time-resolved photoelectron spectroscopy or 

as the gene ralization of the “nanoplasmonic-field micro¬ 


scope” [2l| from surface plasmons to the study of single 
molecules. The key advantage of self-streaking is that 
time-resolved probing of the coherent electron dynamics 
proceeds in the absence of any distorting external NIR 
field. 


The electron dynamics in NeQCeo driven by the UV- 
VIS field are treated within time-dependent density func¬ 


tional theory [22 
elsewhere 26, 27 


[25j . The simulation has been described 
Briefly, we solve the time-dependent 
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Kohn-Sham equations 

idt'ipz{r,t) = I - + Mon + FLit)z + 

J (1) 

where the terms on the right hand side correspond to the 
kinetic energy, the potential of the carbon and neon Is^ 
cores in terms of norm-conserving pseudopotentials [1^ , 
the coupling to the UV-VIS laser field F^it) in dipole 
approximation and length gauge, the Hartree potential, 
and the exchange-correlation (XC) potential for which we 
employ the adiabatic local density approximation 
The valence electron density is n(r,t) = OP- 

Eqs. © are integrated in real space in a cuboid box 
of 29 X 29 X 67 a.u. with a space discretization of 0.4 
a.u. employing a nine-point stencil for the spatial deriva¬ 
tives and a fourth-order Taylor approximation to the 
time-evolution operator. The time-dependent induced 
dipole moment is the first moment of the induced elec¬ 
tron density An(r, t) = n(r, t) — n(r, —oo). We choose a 
UV-VIS field resembling those currently attainable 
[IS El and record the three-dimensional microscopic 
electric field distribution FNF(r, t) evaluated as the nega¬ 
tive gradient of the induced potential V (r, t) — V (r, —oo) 
(Eq. (P)) [Hil. 

The ab-initio near-held distribution serves as input 
for subsequent classical trajectory Monte Carlo (CTMC) 
simulations [s^. For the propagation of the photoemit- 
ted electron the CTMC simulation of streaking is well 
justihed for attosecond XUV pulses with pulse duration 
Txuv short compared to the oscillation period of the 
streaking held, txuv (T^: period of the dipole 

streaking held) when interference effects are negligible 
and for moderately fast electrons (velocity w > 1 a.u.) 
Hi [13. The initial distribution of the photoelectrons 
is given by the neon 2s spatial distribution with an ini¬ 
tial dipolar angular distribution with anisotropy param¬ 
eter /32s = 2. Their initial mean energy is given by 
Ei = huxvY — Ip with 2s binding energy Ip = 48.5 eV 
taking into account spectral broadening by the attosec¬ 
ond XUV pulse (duration txuv = 100 as). The ionization 
probability is assumed proportional to the instantaneous 
XUV intensity. The electrons are propagated in the space 
and time dependent molecular near-held ENF(r, t). Their 
hnal momenta within an acceptance angle of 10 deg from 
the polarization direction are recorded as a function of 
the delay At between XUV and UV-VIS pulses in a 
“streaking spectrogram” (Fig. (T}:). The streaking sig¬ 
nal of each simulation is the delay-dependent shift of the 
center of momentum, AP(Ar). 

Elastic scattering at the cores of the carbon cage 
atoms and inelastic electron-electron scattering along the 
trajectory can by easil y ta ken into account within the 
CTMC simulation [l^ ISSl but do not signihcantly 



delay [fs] 

FIG. 1. (Color online) (a) Laser field (central wavelength 310 
nm, 4 fs pulse length, sin^ envelope, intensity 9x 10^^ W/cm^) 
(b) Dipole moment of NeQCeo- Inset: Near-held (left) and 
near-held of dielectric sphere with the same dipole moment 
(right, Eq. O at t = 6.7 fs). Circles mark the atomic posi¬ 
tions. (c) Simulated streaking spectrogram (taxuv = 100 eV, 
duration 100 as). Black line: center of momentum streak¬ 
ing signal AP(Ar). (d) Enlarged streaking signal after the 
conclusion of the laser pulse (black solid, squares) and time 
shifted dipole moment (red dashed). 


influence the near-held induced shift of the mean mo¬ 
mentum. The most likely source for the uncertainty in 
the simulation comes from neglecting the back-action of 
the core hole on the response of NeQCeo- The time scale 
on which the Cgo spectator electrons will respond to the 
additional Coulomb potential of the ionized neon core can 
be estimated from TD-DFT simulations to be around 10 
a.u. (0.25 fs) [i^, M| starting from the time the electron 
passes the Cgo shell. The ionized electron has then al¬ 
ready reached a distance of ^ 20 a.u. or three times the 
Ceo radius from the shell and has sampled most of the 
(de)acceleration exerted by the molecular near-held. 

During the laser pulse, the dipole moment oscillates in 
phase with the electric held (Fig. [TJ)) while the streaking 
shift AP is in phase with the vector potential (Fig. [1]:), 
i.e. 90 deg shifted relative to the held as in conven¬ 
tional streaking. After the conclusion of the laser pulse, 
the time-dependent dipole moment continues to oscillate 
(Fig. [T]3) with a dominant oscillation period ^ 1 fs corre¬ 
sponding to the transition frequency to the lowest dipole- 
allowed excited states — Eq ^ i eV. The longer beat- 
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iiig period ^ 5 fs corresponds to the energy difference be¬ 
tween different coherently excited states Eg, — ife ^ 1 eV. 
This time dependence of the dipole moment indicates 
changes in the induced density as well as in the alter¬ 
nating dipolar near-field (Fig. [TJ) inset). The nano-scale 
electric field distribution resembles the sum of the laser 
Fz,(<) = FL{t)z and of the field of a dielectric sphere with 
dipole moment d{t) = d(t)z, 


FNF(r, t) 


FL{t) + 


-d(t)/R^ for 

—d(t)/r^-I-for 


r < R 


r > R 

( 2 ) 


where R = 9 a.u. matches the field enhancement and 
near-held of the effective electronic surface somewhat 
outside the ionic positions at 6.66 a.u. (Fig. [T}d inset) 
For time delays exceeding the duration of the UV- 
VIS pulse, the streaking sprectrogram encodes the in- 
huence of the time-dependent dipolar near-Held on the 
emitted electron. We Hnd that the streaking modulation 
AP(Ar) traces the oscillations of the dipole moment to 
a surprisingly good degree of approximation. The molec¬ 
ular dipole moment d{t) can be reconstructed from the 
streaking signal AP(Ar) as 


d{AT) = AAP{AT + At) (3) 


with a scaling factor A ~ 190 a.u. and a time shift At ~ 
0.36 fs (Fig. [TJ:). 

The mapping of the molecular dipole moment onto the 
streaking signal (Eq. [3]) is remarkably robust under pa¬ 
rameter variation. We numerically obtain a linear de¬ 
pendence of the streaking signal on the pump laser Held 
amplitude while its temporal shape remains almost un¬ 
changed when varying the intensity over 2 orders of mag¬ 
nitude (Fig.[3Ji). Similarly, the streaking time signal does 
not depend sensitively on details of the near-Held. When 
either changing the molecular orientation, or removing 
the central neon atom, or even replacing the Held dis¬ 
tribution from TD-DFT by the polarized sphere Eq. [2] 
(Eig. Od), the mapping of the time-dependent induced 
dipole moment to the streaking signal closely agrees with 
that of the reference calculation (Eig.[TJ:). Averaging over 
molecular orientations and, thus, over the anisotropy of 
the molecular polarizability, we Hnd that the orientation- 
induced dephasing is negligible on the time scale of our 
simulation. 

We now analyze the properties underlying the remark¬ 
ably stable mapping of the molecular dipole moment onto 
the streaking signal AP, employing the simpliHed Held 
distribution (Eq. [2]) for convenience. The Held distribu¬ 
tion in the absence of the laser pulse (Eq. [3] with laser 
FL{t) = 0 and d{t) > 0) is symmetric under reflec¬ 
tion (z —i> —z) at a plane perpendicular to the dipole 
polarization. Correspondingly, the electrostatic poten¬ 
tial is antisymmetric and is zero both at the origin and 
at inflnity. An electron emitted from the center of the 
Ceo will emerge at large distances with its initial energy 



FIG. 2. (Color online) (a) Intensity dependence of streaking 
signal; reference calculation from Fig. [TJ; (9 x 10^^ W/cm^, 
black solid); lower intensity (1 x lO^^W/cm^, red dashed, 
x3), higher intensity (1 x 10^®W/cm^, blue dotted, xO.3) 
each scaled by streaking field amplitude. Streaking momen¬ 
tum shift is 0.007 a.u. (0.17 eV), 0.19 a.u. (0.51 eV), 0.67 
a.u. (1.82 eV) respectively, (b) Dependence on near-held; ref¬ 
erence calculation from Fig. [TJ: (black solid), different orienta¬ 
tion of NeQCeo (green dashed) and of C6o(blue dash-dotted), 
sphere model (Eq. [2] turquoise dotted), (c) Decay of dipole 
moment; reference calculation from Fig. [TJ: (black solid) for 
one orientation, orientation averaged (green dotted), orienta¬ 
tion average with exponential decay with decoherence time 
tc = 20fs (details see text, red dashed). 


as the deceleration inside the sphere is exactly canceled 
by the acceleration in the enhanced near-held outside for 
a static Held distribution (Eq. [2|). The observed delay 
time dependent momentum shift AP and, thus, the self¬ 
streaking relies on non-adiabatic corrections to the static 
Held. The time scale for the traversal of the molecular 
near-Held extending over ~ 2R for a photoelectron with 
velocity vq — Y^2(lia;xuv “ Ip) ~ 2 a.u. is comparable to 
a fraction of the dipole oscillation period ~ 7^/4, over 
which the Held changes from zero to its maximum value, 
2R/vo ~ Td/4:. This qualitative argument can be made 
quantitative by further simplifying Eq. [2j Accordingly, 
the streaking momentum shift for an electron ionized at 
birth time At along the laser polarization direction z is 
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approximately given by 


AP(Ar) = / Fjsif {z{t),t) dt 
J Ar 

doTd . f R \ . 

R^tt V "^oTd J 


sin TT 


R + 2woAr \ 
voTd ) 


+ 


do i ( R + voAt\ 

sin + 


2Tr^R'^ 


VoTd 

iS) i^’^) 


• (4) 


Here we have assumed for simplicity (i) the dipole mo¬ 
ment oscillates with a single Fourier component d{t) = 
do sin(27rt/Td) and (ii) a straight-line trajectory z(t) = 
vo{t — Ar) of the emitted electron. The first term rep¬ 
resents the motion inside the cage and the second term 
the motion in the enhanced near-held. Eq. 2] effectively 
predicts values for A and At (Eq. [31) in good agreement 
with the simulation (see Fig. |31 d). The largest momen¬ 
tum shift is incurred for electrons emitted shortly after 
the extrema of the dipole moment (Fig. |3l times ti and 
to). These trajectories are accelerated inside the cage and 
reach the molecular surface at the time when the dipole 
moment changes its sign so that they are further accel¬ 
erated by the enhanced near-held. For electrons born at 
a time when the magnitude of the instantaneous dipole 
moment increases (Fig. El ^ 2 ), the self-streaking held is 
stronger but the held inside and outside the cage have 
different signs and therefore partially compensate each 
other leading to a reduced momentum shift. The present 
self-streaking protocol is independent of how the dipole 
excitation is initiated. Alternatively to the moderately 
intense VIS-UV pulses employed here, lower frequency 
pulses of higher intensity could be used as well [4^. 

Eventually, the coherent electronic dipole excitation 
will dephase by coupling to vibrational degrees and decay 
by electron-electron scattering or radiative relaxation, 
the latter of which occurring on the nanosecond scale. 
For Ceo, the fastest dephasing and dissipative process is 
estimated to be coupling to the lattice degrees of free¬ 
dom. An estimate for its time scale is given by the os¬ 
cillation period of the C=C stretching mode of ^ 23 fs. 
This agrees well with estimates of electron-phonon cou¬ 
pling from the experimental spectra and theory 4^47|. 
Coupling to other degrees of freedom will lead to damp¬ 
ing of the dipole oscillation amplitude and, in turn, of 
the streaking amplitude (Fig. |2t). Self-streaking, thus, 
affords the opportunity to study dephasing and dissipa¬ 
tive processes directly in the time domain and benchmark 
emerging ab-inito descriptions of non-adiabatic electron¬ 



time [fs] 

FIG. 3. (Color online) Enlarged time interval [5.4 fs — 6.9 fs] 
from Fig. [2] (a) Dipole moment; reference TD-DFT cal¬ 

culation from Fig. [T] (black solid), single-harmonic fit (blue 
dotted), (b) Streaking signal; reference TD-DFT calculation 
from Fig. [T] (black solid), simplified sphere field distribution 
(Eq. [^ green dashed), analytical model (Eq. [4] blue dotted), 
(c) Force along electron trajectories starting at ti = 5.56 fs, 
t 2 — 5.74 fs, to = 6.02 fs (sphere field distribution Eq.[2j. (d) 
Momentum shift accumulated along these trajectories. 


ion coupling [4^ and electron-electron scattering pro¬ 
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cesses 

In conclusion, we have demonstrated that dipolar ex¬ 
citations in the Ne^Ceo molecule can be investigated 
with attosecond resolution employing a novel variant 
of the attosecond streaking protocol involving intrinsic 
rather than external oscillating strong fields. This “self¬ 
streaking” maps out the real-time modulations and de¬ 
cay of the time-dependent molecular near-field and dipole 
moment. For endohedral fullerenes, photoionization of 
the central rare-gas atom provides a localized source of 
distinguishable fast streaking electrons. Similarly, for 
studying processes in more complex systems such as or¬ 
ganic molecules containing several atomic species, contri¬ 
butions from different locations within the molecule can 
be disentangled by the final electron energy. We may en¬ 
vision time-domain visualization of charge-transfer pro¬ 
cesses [13, [EH and hot electron dynamics in small clusters 
and nanoparticles 52 -5^ with sub-femtosecond resolu¬ 
tion. 
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